library(scRNAseq)
library(velociraptor)
library(scuttle)
library(scran)

sce <- HermannSpermatogenesisData()
sce <- sce[, 1:500]

sce

sce <- logNormCounts(sce, assay.type = 1)

top.hvgs <- sce |>
  modelGeneVar() |>
  getTopHVGs(n = 2000)

velo.out <- sce |>
  scvelo(subset.row = top.hvgs, assay.X = 'spliced')
